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Introduction 


■  High  pressure  heterogeneous  combustion 
systems  typically  operate  trans-  or 
supercritical 

□  Liquid  rocket  engines 

□  Gas  turbine  combustors 

□  Diesel  engines 

■  Quantitative  assessments  or  validated 
theories  characterizing  propellant  injection 
near  and  supercritical  conditions  are 
currently  limited 


SSME 
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Introduction 

■  Because  these  systems  operate  at  near  and  supercritical 
conditions,  it  is  difficult  to  model  the  injection,  mixing, 
atomization,  and  vaporization  using  traditional  CFD 

System  becomes  fully  3-D 

Material  and  transport  properties  change  significantly 
□  Moves  from  2-phase  to  1 -phase 


□  Complex  mixture  effects 
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Introduction 

■  This  research  attempts  to  address  these  limitations  through  the  use  of 
MD  as  an  alternate  simulation  tool 


■  Diatomic  liquid  nitrogen  (N2)  is  simulated  as  a  non-reacting  substitute 
to  other  diatomic  cryogenic  liquid  propellants/oxidizers,  i.e.  H2  or  02 

Experimental  database  exists 


■  The  simulation  consists  of  two  cylinders 

Smaller  diameter  cylinder  representing  one  injector  through  which 
simulated  liquid  nitrogen  flows 

□  Larger  diameter  cylinder  of  gaseous  nitrogen  into  which  the  simulated 


i  c 


Injector  Element 


“Combustion  Chamber” 
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Introduction 

■  Advantages  of  the  MD  model: 

Based  purely  on  first  principles  ( F=ma ) 

Can  theoretically  span  all  conditions;  no  need  to  track  phase 
boundaries 

Geometric  assumptions  do  not  need  to  be  made 

Equation  of  state,  material  and  transport  properties  are  not 
needed,  and  can  actually  be  calculated  outputs 

□  Mixture  effects  automatically  incorporated 
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Fundamental  Concepts  of  MD 


■  Dynamics  of  interest  is  contained  within  solution 
to  the  classical  N-body  problem 

■  MD  is  a  numerical  tool  used  to  solve  Newton’s 
Laws  of  Motion,  applied  to  motions  of  individual 
atoms  or  molecules 

□  F  =  ma  is  the  only  equation  to  solve 

□  Forces  on  atoms  derived  from  interatomic  potentials 
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Fundamental  Concepts  of  MD 

■  Lennard-Jones  12-6  Potential  For  Nitrogen  Atoms: 


■ 
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Fundamental  Concepts  of  MD 


Equations  of  Motion 

N0 


m 


a 


=4=1  /; 


j=i 


For  Lennard-Jones  Potential 
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Solve  Numerically  Using 
Finite  Difference  Velocity- 
Verlet  Algorithm 


r(t  +  At)  =  r(t )  +  A  tv(t)  + 


A  t2a(t) 


v(t  +  At  /  2)  —  v  ( /  )  + 


Ata(t) 


v(t  +  Ar)  —  ■+■  At  /  2)  + 


A  tu(t  +  Ar) 
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Fundamental  Concepts  of  MD 

■  Addition  of  Bond  Constraint  to  Simulate  Diatomic  Nitrogen  Molecule 


rab  =  rab  •  rab  =  constant 


dt 


=  2rab  ■  Vab  =  0 


RATTLE  Algorithm: 


r(t  +  At)  -r'(t  +  Ar)  + 
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molecule  simulated  separately 
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Simulation  Methodology 

■  Computational  Domain  and  Processor  Setup: 


□  Each  processor  assigned  to  handle  particles  within  a  specific 
geometric  region 

As  particles  move,  they  are  transferred  between  processors 
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Periodic  Boundary  Conditions: 

Flow  Direction 


Cell  Boundaries  Each  2.5a  in  Length 


10 


International  Symposium  on  Energy  Conversion  Fundamentals,  Istanbul,  Turkey,  June  21-25 ,  2004. 


Wall  Model 

■  Constant  Temperature  Wall 

■  Diffusely  Reflecting 

□  Molecules  reemitted  from  wall  in  random  direction 
with  random  velocity  based  on  Maxwellian  distribution 
scaled  to  the  wall  temperature. 

Removes  heat  generated  by  viscous  dissipation 
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Simulation  Methodology 
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Initialized  Liquid 
Nitrogen  Flow 


Initialized  Gas 
Chamber 
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Results 


Radial  Velocity  Profile 
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Results 
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Demonstrating  Rayleigh  Breakup 
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Results 
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Results 


124  K,  ~3.0  MPa,  Re=4.04,  We,=41.53,  We  =7.06 


100  ps 


125  ps  175  ps  200  ps 


Source:  Mayer,  W.  O.,  and  Telaar,  J.,  “Investigation  of  Breakup  of  Turbulent 
Cryogenic  Variable-Density  Jets.”  Atomization  and  Sprays,  12,  651-666,  2002 
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I  Demonstrating  Second  Wind-Induced  Breakup 
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Results 
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Ambient  Gas  Pressure  (MPa) 

■  A  decrease  in  breakup  length  with  an  increase  in 
ambient  pressure  is  observed 

■  The  difference  in  breakup  length  between  the 
higher  ambient  pressure  cases  is  also  less, 
indicative  of  higher  pressures  having  less  of  an 
effect  on  breakup  length 
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Animation 

83  K 

~0.2  MPa 

Re=(a)0.87,  (b)1.70 
We|=(a)0.50,  (b)1 .90 
Weg=(a)0.005,  (b)0.021 

uavg=(a)33.8  m/s 


(b)66.1  m/s 


L/D=60 


Total  Simulation  Time: 
1.2625  ns 


(a)  Case1new.avi,  (b)  fgrav2xspray.avi 
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Supercritical  -  Comparison  to  Raman  measurement 


■  4  MPa,  123  K 


*Branam,  R.  and  Mayer,  W.,  “Characterization  of  Cryogenic  Injection  at  Supercritical  Pressure,”  J.  Prop,  and  Pow.,£foA.  19, 
No.  3,  May-June  2003,  pp.  342-355. 
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Supercritical  -  Computed  radial  density  profiles 
■  Target  conditions,  6  MPa,  123  K,  40  m/s 
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Supercritical  -  Evolution  of  computed  density  profiles 

■  Developing  Flow 

□  At  a  given  time 

□  Average  of  five  independent  simulations 

■  Supercritical  conditions,  4  MPa,  123  K 
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Summary 


■  Subcritical  Injection  Simulations: 


□  Rayleigh  breakup  is  reproduced,  resulting  in 
formation  of  droplets  due  to  capillary  instabilities  in 
the  jet  driven  by  surface  tension 


□  The  onset  of  aerodynamic  effects,  including  the 
second  wind-induced  breakup  regime,  can  also  be 
seen  in  the  cases  of  higher  gas  pressures 
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Summary  (II) 


■  Supercritical  Injection  Simulations: 


□  MD  simulations  match  Raman  measurements 
near  injector  exit  (potential  core) 

□  Simulated  flow  still  transient  further 
downstream 


24 


International  Symposium  on  Energy  Conversion  Fundamentals,  Istanbul,  Turkey,  June  21-25,  2004. 


Future  work 


■  Resolve  wall  B.C.  rotational  accommodation 

■  Increase  size  of  simulations 

Generating  injected  particles 

□  Increase  tube  diameter 

□  Increase  Re,  decrease  Kn 

■  Future  work  to  focus  more  on  ethanol  and  other 
hydrocarbons 

■  Extend  geometry  to  include  coaxial  tube 
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